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Abstract. We investigate the effects of hyperbolic hydrodynamical mixing on the 
reaction kinetics of autocatalytic systems. Exact results are derived for the two 
dimensional open baker map as an underlying mixing dynamics for a two-component 
autocatalytic system, A + B — > 2B. We prove that the hyperboliticity exponentially 
enhances the productivity of the reaction which is due the fact that the reaction 
kinetics is catalyzed by the fractal unstable manifold of the chaotic set of the reaction- 
free dynamics. The results are compared with phenomenological theories of active 
advection. 
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1. Introduction 



The advection of chemically or biologically active particles in imperfectly mixing chaotic 
hydrodynamical flows received considerable recent interest |], ||, ||, |, [5], || |7], |8|, By 
chaotic we understand a smoothly time-dependent, in the simplest case periodic, flow in 
which the advection dynamics of tracers is chaotic. These tracers undergo chemical or 
biological interactions while being advected by two-dimensional flows. The distribution 
of the products along a filamental fractal seems to be a general feature of such reactions. 

Here we consider tracers as point-like particles which are assumed not to have 
any influence on the fluid flow itself. The reaction model is of kinetic type where the 
reactions take place upon particle-particle 'collisions': when particles of type A and B 
are closer to each other than a reaction range a, then particle A changes to B according 
to the autocatalytic scheme A + B — > 2B. The reaction model is discrete in time, i.e., 
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a reaction can only take place at integer multiples of a reaction time r. In intervals 
between successive reactions the particles are passively advected by the flow. 

The passive advection dynamics has two essential parameters that influence the 
reaction kinetics of the process: a decay or contraction parameter A and the fractal 
dimension Do of the set around which products accumulate. In open flows this set 
turns out to be the unstable manifold of a nonescaping chaotic set characterizing the 
passive advection [Q, ||], and the role of the contraction parameter is played by the 
average negative Lyapunov exponent of the Lagrangian dynamics which expresses the 
convergence towards the chaotic set. Due to incompressibility, the positive and negative 
Lyapunov exponents have the same absolute value. Particles are advected away from 
any neighborhood of the chaotic set along the unstable manifold which leads to an 
exponential decay of the passive particles with an escape rate k. This is, however, not 
a new dynamical parameter since it is related to the previous ones via Dq »s 2 — k/X, 
flO| , |H| . Here we concentrate on slow reactions whose reaction velocity v r = ojr is 
much slower than the characteristic velocity of the passive advection dynamics. 

In this section we present a heuristic argument for the total area of the product 
particles as a function of time. Since the basic ingredients of the argument are 
the existence of a fractal backbone and the contracting dynamics acting towards the 
backbone, we do not need to actually specify the two-dimensional hydrodynamical flow 
itself just to postulate the fractal backbone and the contraction property. We consider 
the simplest filamental fractal which is the direct product of a Cantor set of dimension 
(Dq — 1) and a line segment of length L so that the total dimension is -Do- In addition, 
a contraction of rate A is postulated towards the fractal filaments, see Figure |I[ This 
simple theory can also be considered to hold for reactions taking place along spatially 
fixed fractal catalysts when there is transport of products, resulting in an effective 
contraction. 
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Figure 1. Schematic drawing of the emptying and fattening process surrounding a 
single filament. A line segment of length L ( a)) is covered by material B ( b)). The 
width of this stripe is first shrunk ( c)) due to outflow, then it is fattened by a sudden 
reaction ( d)). 



At time n, right after a reaction takes place, filaments are assumed to be covered by 
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material B in stripes of average width eW. The complement of these strips is assumed 
to be covered by material A everywhere. During the next period of length r no reaction 
takes place, only contraction, which reduces the widths e^ 1 ' by a factor of exp (— At). 
The effect of the reaction is simply an 'infection' by B particles of the A neighbors in 
circles of radii a which in case of a filamental geometry leads to the broadening of the 
stripe width by ga, g = 2. Thus we obtain a recursion relation expressing the stripe 
widths ght after the (n + l)-st reaction in terms of e^, right after the n-th 

reaction: 



e (n+i) = e -^ £ (n) + g(J ^ 
A7 (2) 



This simple recursive dynamics possesses a single fixed point attractor 

9° 



l-e~ 

which is reached for n — > oo. This shows that the product distribution reaches 
saturation: an active steady state sets in after some time in which the activity exactly 
cancels the material loss due to contraction. 

Although there are infinitely many filaments in the fractal backbone, if they are 
covered with stripes of non-zero width, the number of such covered filaments becomes 
finite due to the overlaps among them. It follows from the properties of fractals that 
the number N(e) of filaments covered by width e is 

/ e \ -(A)-i) 

N(e) ~ ( z j . (3) 

The total area A of these covered filaments is then 

A(e) = eLN(e) ~ e 2 - D °L D °. (4) 

Next we apply this result with e = at time n and with e = e^ n+1 ^ at time n + 1. Note 
that we are then working on the egde of the validity of fractal scaling since below the 
scale of or the product distribution is two-dimensional. The consequences of 

the crudeness of this approximation will be analyzed later, in view of the exact results 
obtained for the baker map. Using one finds a recursion for the area Ab of B 
particles. It follows from ([[]) that 



A (n+1) = I e -X T 



A (n) 



l/(2-D ) 



2-D 



+ gaL 2 ~ D o \ . (5) 



In the case of exactly parallel filaments the factor g is identically g — 2 before the 
stripes start to overlap. After overlap takes place, which is unavoidable on sufficiently 
fine scales, the value of g may however differ from 2. In this approach, we assume that 
the overlap effects can be taken care by introducing an average geometrical factor g, 
which is independent of the reaction parameters a and r. If the proportionality factor 
in (|4]), which is the Hausdorff measure of the fractal set, is different from 1, it may 
also be incorporated in the geometrical factor g. The form of recursion (|5|) is a lot 
more complicated than recursion ([!]) for the stripe widths due to the fractal character 
of the backbone. It shows that the 'microscopic' width dynamics and the 'macroscopic' 
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area dynamics are qualitatively different since the latter is nonlinear. Nevertheless, 
possesses a fixed point 



ga/L 



2-Do 



which is the only attractor of the chemical dynamics. 

It is worth taking a continuum limit in which the reaction time r and distance a 
both go to zero, while the reaction velocity v r = cr/r remains finite. In this limit a 
differential equation is obtained for the B area As{t) as a function of the continuous 
time t in the form 

A B = -A(2 - D )A B + g(2 - D )v r L*$> {A B )~ P (7) 
where the exponent f3 is a unique positive expression of the fractal dimension 

"TO < 8 > 

The singular-looking reaction term on the right hand side (rhs) of (|?p expresses the fact 
that the smaller the area, the larger the production rate becomes, just as one expects 
for reactions on a fractal catalyst. 

Both the discrete-time and the continuous-time area dynamics possess a fixed point 
steady state. The corresponding product area of the latter is 

Note that this is much larger than the product area along a single non-fractal line 
(obtained as the special case Dq = 1) since the expression in the paranthesis is much 
less than unity in view of our assumption of the slowness of the active dynamics. This 
clearly indicates that the presence of a fractal catalyst makes the reactions much more 
efficient. The important caveat is that the fractal backbone is not introduced 'by hand' 
as an imposed mixing rule, it is rather generated naturally by the chaotic advection 
dynamics of the fluid flow transporting the tracers. The chaotic set is the union of 
unstable bounded orbits. The reactions take place along the outflow curves from these 
orbits (the unstable manifold). For a time periodic flow the unstable manifold moves 
periodically in the plane of the flow (but does not change its fractal dimension). This 
is, however, not a basic difference since on a stroboscopic map the manifold appears to 
be still. 

In order to illustrate and also check the above general heuristic arguments, in 
this paper we consider two variants of the baker map, which describe a discrete mixing 
dynamics with hyperbolic chaotic properties. These maps are simple enough to facilitate 
an exact analytic approach and yet they incorporate the generic features of mixing from 
realistic flows. Within this model we are able to derive exact expressions for the area 
dynamics and thus to identify the geometric factor, which was not possible within the 
phenomenological theories of active chaos. 

Our paper is organized as follows: in Section ||] we introduce a simple variant of 
the baker map on which we superimpose the chemical activity, introducing thus a novel 
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reactive dynamical system, the Active baker map. Then we derive the exact stripe 
dynamics for this map, including the reaction equation describing the area dynamics of 
the reactants in the system. In Section |^ we compare the exact expressions with the 
heuristic results given in the Introduction, identifying the exact effects of the fractal 
geometry on the geometric factor. In Section || we consider a second variant for the 
baker map, and show the robustness of the heuristic results. Section |^ is devoted to 
discussions and outlook. 

2. The active baker map 

The passive (reaction-free) baker map is probably the simplest dynamical system which 
exhibits chaotic behaviour of purely hyperbolic type. Its simplicity allows for analytic 
manipulations, and thus it plays from a theoretical standpoint a crucial role, comparable 
to the Ising model from statistical mechanics. There are several variants to this map, 
and we shall consider two of them in more detail. Let us give first one variant of the 
baker map and then define the active version based on this, then in Sec. ^ introduce 
another variant with the associated active counterpart, and then finally compare the 
results based on the actual realizations of the maps. 

Assuming that we start from a completely filled unit square, one step of the baker 
map consists of the following: 1) the lower (upper) half of the unit square is compressed 
in the x direction by a factor of a < 1/2 around the fixed point in (0, 0) ((1, 1)), and then 
2) the compressed stripes are stretched along the y direction (around the same fixed 
points) by a factor of b > 2. For b > 2 there will be 'material' reaching over the unit 
square, which is clipped and discarded. In this case the map is modeling an open flow. 
By taking the unit square as the phase space of the baker map, the characteristic length 
L of the dynamics is chosen to be unity. When ab — 1, the baker map is area preserving 
in the sense that its Jacobian is uniformly unity, however, it is an open chaotic map 
for any b > 2, the chaos being of transient type [pLlJl . Mathematically, the image {x',y') 
through the baker map of a point (x, y) from the unit square is given by: 

x' = ax + (l-a)6(y-l/2) , x E [0, 1] 

y> = by-(b-l)0(y-l/2) , ye [0,1] 1 } 

where 8(x) is the Heaviside step function. The map is not defined outside the unit 
square (i.e. L — 1) and thus chemical reactions can only occur inside the unit square. 
In the following we restrict ourselves to open area preserving systems, i.e., to the regime 
b = 1/a > 2. Repeated application of the map will result in a sequence of vertical stripes 
of various widths which form the hierarchy of a one dimensional Cantor set in the unit 
interval along the x-axis. Contraction towards the vertical filaments is governed by the 
Lyapunov exponent 

A = -lna>0. (11) 

Since due to the construction, in the vertical direction there is no structure, the map is 
effectively one dimensional. The fractal dimension Dq — 1 of the Cantor set on which 
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the stripes reside, is simply calculated to be: 

In 2 

Dn-1- 



lna 



(12) 



The decay of the area covered by the material left in the unit square is exponential, 
with a rate given by the so called escape rate, k: 

K = -m(2a), (13) 

so that D = 2 — k/\ holds for the fractal dimension of the stripes in the unit square. 
Initially, the empty space between the stripes is imagined to be filled with 'material' A, 
while the stripes contain material B. 

We call the baker map active when the two materials react with each other. We 
consider autocatalytic reactions A + B — > 2B (the same as in the Introduction), with the 
reactions happening along the (vertical) contact lines. This reaction can alternatively 
be considered simply as the surface-growth of material B, and leads to a broadening of 
the stripes relative to those of the passive baker map. 

We assume that the reactions happen instantaneously and simultaneously in the 
system, after every r steps of the reaction-free baker map. Time is measured in units 
of the baker-steps. Thus, for r = 1, the B surface grows after each baker-step by an 
amount of a on each side of a stripe. For r = 2, two reaction-free baker-steps are taken 
before one instantaneous reaction occurs, etc. 
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Figure 2. Two consecutive steps of the baker map and two reactions (r = 1). The 
bands of width a become occupied by B in each reaction. The material hanging over 
the unit square is discarded. 

For the simplicity of writing, let us introduce the following convention: passive steps 
will be called baker-steps, and the combination of r baker-steps followed immediately 
by one reaction will be called one active-step. 

If the gap between two stripes is shorter than 2cr, the gap disappears during reaction, 
and the stripes merge into a single one. In the absence of reactions (a = 0) the steady 
state (defined as the asymptotic state t — > oo) is characterized by a zero Lebesque- 
measure distribution of material B, covering the unstable manifold of teh chaotic set of 
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dimension D of the baker map. In the case of the active baker map, however, there 
are new I?-areas of width a appearing after periods r in the system, and thus a sort of 
dynamical balance is expected to set in between the escape dynamics and the growth 
dynamics. 

2.1. Stripe hierarchy in a simple case: one baker-step per active-step 

In spite of the fact that the baker map is one of the simplest maps that exhibits chaotic 
dynamics, when chemical activity is also introduced, the reaction kinetics of the system 
becomes rather involved. 

Let us first start with the simpler case of one reaction per baker-step, i.e., r = 1. 
Figure ^| shows the succession of two such active-steps. The black bands show the B 
material coming from the first reaction, while the dotted light gray bands are the B 
material emerging from the second reaction. From the point of view of the reaction 
kinetics, the distribution and the amount of material A (the gaps) is crucial. If an A 
stripe is less wide than 2a then it will fill up in the next reaction and disappears. During 
the baker-steps the hierarchy of gaps builds up with thinner and thinner gaps appearing. 
From a certain level on, the reactions are putting a stop to the hierarchy growth. The 
question is how in the long time limit the picture stabilizes and the system settles into 
a stationary state. Every gap is bounded from left and right by material B. Since the 
contraction by a is a point transformation (around the two fixed points at the edges), 
each gap will shrink by a factor of a and it will also get closer to the corresponding 
fixed point by the same factor a. After the reaction the gaps will lose 2a from their 
width. After one (n = 1) active-step we get one gap of width 1 — 2a — 2a; after two 
active-steps we obtain three gaps: one of width 1 — 2a — 2a, and two gaps, each of width 
a(l — 2a — 2a) — 2a; for n = 3, we have a total of 7 gaps, three just like in the n = 2 
case and 4 new members each of width a[a(l — 2a — 2a) — 2a] — 2a, etc. An important 
observation is that if gaps appear at a certain level of the hierarchy, there always will 
be gaps of precisely the same width and position in the unit square at any later time: 
during the next active-step, a particular gap will map precisely down into the location 
of two other gaps of the given level. The largest gap (the top level of the hierarchy) is 
always located in between the points xi = a + a and x r = 1 — a — a. During the next 
active-step this particular gap is mapped into two lower level gaps, however its original 
position and width is taken over by the new gap appearing as the images of the edges 
at x = and x — 1 of the unit square, and modified by a due to the reaction. 

Thus the hierarchy being formed can be described as follows. After n active-steps, 
there will be gaps of type n (i.e., the hierarchy will have n levels). If Wi denotes the 
width of a gap on level j, and denotes their number (or multiplicity), then one can 
write: 

w? = (1 - 2aV'- 1 - 2a ] —^- , m[ j) = 2 3 '" 1 . (14) 

1 — a 

Observe that the first term on the rhs of the first equation is exactly the width of the 
gaps in the reaction free baker map [a = 0). The second represents the total amount of 
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area-decrease of each gap in j reaction steps: — 2cr(l + a + a 2 + ...a^ 1 ). The total area 
occupied by the gaps (or material A) after n active-steps is given by 

Af = ± m? W ? = 1 - (2a)" - ( 2 » - 1 - al^f} (15) 

The area occupied by material B is simply the complementer of the above in the unit 
square: 

A f = 1 - A f = (2a)" + ( 2 » - 1 - ai^^) • (16) 

When n —>■ oo, .4.^ becomes negative, and later both areas formally diverge which 
indicates that this process has to stop after a finite number of steps. This happens by 
the closing of the smallest gaps due to the reactions in the n*-th step: no further gaps 
can be created from a certain level down in the hierarchy because the reactions are of 
non-zero, finite width a. By definition, the thinnest gaps of an n-level hierarchy are the 
gaps at the lowest (n-th) level of width w^K The condition for the lowest level to be 
filled in n* active-steps is w^ 1 ^ > and w[ n < 0. From this, and Eq. [14], we obtain 

* _ f [u] + 1 , if u is non-integer , 
\ u , if u is an integer , 

with 

l n (l + (l-«)(l-2a) \ 

« = ; , (18) 

a 

([x] denotes the integer part of x). For n > n*, < and there are no more gaps 
at and below this level. Thus, even if the active-steps are repeated indefinitely, the 
hierarchy of gaps stops at j = n*. Since the gaps at higher level map into gaps at 
lower level, a steady state is reached after n* steps, and the whole material distribution 
becomes stationary. The fact that a stationary state is reached as a balance between 
the emptying dynamics and the reactions was numerically observed in all the previously 
studied dynamical systems (e.g., flow in the wake of a cylindrical obstacle) |2|, U, 
corroborating the generic character of the stationary state. 



2.2. Stripe-hierarchy for multiple baker-steps per active-step 

In this case there are r baker-steps taken before an instantaneous reaction occurs. Let 
us describe the hierarchy of gaps built up during n active-steps. After the first r baker- 
steps we obtain the usual reaction-free baker hierarchy with r different types of gaps (a 
total of 2 T — 1 gaps). After the first reaction there will be therefore = 2 k ~ l gaps of 
type k whose width is = a fc_1 (l — 2a) — 2o with k = 1,2, .., r. Following the next 
active-step, one obtains 2r types of gaps. The first r types of gaps are generated by the 
k times iterated images of the endpoints x = and x — 1 (A; = 1, 2, .., r) and modified 
by the reaction. These have multiplicities widths and positions exactly as those after 
the first active-step. The second r types of gaps are r times iterated images of the gaps 
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coming from the first active-step. The widths and multiplicities of this second set are 
given by = a T w^ — 2a, and = 2 r+fc ~ 1 on each level k — 1, ..,r. In general, 
if we perform n active-steps, there will be a total of 2 nr — 1 gaps generated with the 
following widths and multiplicities: 



W W = a k -\l - 2a) - 2a , mft = 2^ 



w 



W = a T w ( i~ 1) - 2a , mf = 2^>' 



-i 



(19) 



where k = 1, ..,r and j = 2, ..,n. Analogously to the r = 1 case, during the iteration, 
the gap-hierarchy changes by the addition of finer gaps, and the previously created gap- 
structure is not modified by the new reaction. However, in contrast to the r = 1 case 
the new gaps resulted from the new reaction form a sub-hierarchy of r levels. In the 
expressions ([H|) j denotes the main levels of the hierarchy (related to the state just 
after a reaction) while k quantifies the hierarchy on the subtree level (created by the 
in-between- reactions baker-steps). Recursion flT9]) is solved easily, and the result is: 

j = l,...,n. (20) 



= (l-2a)a^- 1 ^ +k - 1 ~2a^ 



rn 



k 

Obviously the widths are decreasing functions of both j and k (down the hierarchy) as 
also seen from (^0|). Therefore, after n active-steps the thinnest gaps will have the width 
of ?4 n) = (l-2a)a nT - 1 -2a(l-a nT )/(l-a T ) and multiplicity of = 2 nT -\ Another 
important consequence of the monotonic behaviour is that the widths found within a 
level j are all smaller than the widths found within the level j — 1. This statement is 
easily proven by comparing the thinnest and thickest gaps of the neighboring levels. 
Thus, the total area of material A after n active-steps is computed as: 

^ = tt^M J) = 1 - (2a)- - (2- - 1 - a T (2- - 1) \~_ { ^y ) 7^7 ■ ( 21 ) 

For r = 1 we indeed recover Eq. fll5|) . The total amount of material B left in the system 
will again be the complementary of ( pl|) in the unit square, i.e., 

Af = (2ar + (2- - 1 - a' (2T - 1) \~_ ( ^r ) 7^7 ■ ( 22 ) 

As n (time) increases, there will be a critical n* = n*{r) value at which gaps start to 
disappear for the first time. Of course, Eq. (|22|) is valid only in the regime where no 
gaps have been filled yet. This means that for n = n* — 1 all the gaps of the hierarchy 
(POP exist, but for n = n* some of the gaps have disappeared. Let us calculate this value 
and see how a stationary state is reached. Since the thinnest gaps of width must 
have disappeared at n = n* (i.e., < 0), but they were present at n = n* — 1 (i.e., 

w (ji*-i) > 0), n* is calculated as 

n * _ j [«] + 1 , if u is non-integer , 
[ u , if u is an integer , 



where 



l n (l + ft-«0(l-2a) 
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This gives the start of the filling in terms of the number of active-steps. However, it is 
not necessary that only the thinnest gaps of n* have been filled. In principle, all the 
gaps of level n* can be filled because, as observed above, the widths of gaps at a higher 
level are all larger than those of a lower level. How many gaps of level n* are filled does 
not follow from (p4|). 

Performing another active-step, we may still obtain a partially filled sub-hierarchy 
on the main level n* + 1, and similarly in the next step, etc. However, since the widths 
of the gap additions is strongly decreasing, this process will eventually stop. Once 
stopped, the whole material distribution becomes stationary and the system is settled 
into a steady state. 

The system reaches the steady state when also the thickest {k = 1) gaps are filled 
on the next main level. This happens at m* for which Wi > if n < m* , but ?«{ m < 0. 
This leads to: 

\ [s] + 1 , if s is non-integer , 
m = < . (25) 

Is, it s is an integer , 



where 



, n 1 (l-q^)(l-2tt) 
m 1 ~ a T 2a 



(26) 

Next we compare m* and n* in order to see how many steps are needed to reach 
the steady state, once a filling begins at n* . From Eqs. ( f24|) and fl26|) it is obvious that 
s > u (the numerator on the rhs of fl26|) is larger than that of fl2"4|) as r > 1). Let us 
introduce the notation: 



(1 - a T )(l - 2a) 

v = . 

2a 

The quantities s and u then become: 



(27) 



1 + ^^, and u= 1 - + l ^±^. {21 
Th(i) r r in(i) 



Simple manipulations lead to: 




a+v 



„=1--[1 + _^.| <i (29) 

This means that we have the bounds u < s < u + 1, when r > 2. (For r = 1, s and u 
are identical.) This implies that 

m* e {n*,n* + 1} , (30) 

i.e., the steady state is reached either on level n*, or on the next level, n* + 1. 

When m* = n*, Eq. fl2~T|) correctly describes the time evolution of material A over 
the whole range n = 1,2, ..,n* — 1, since there are no partially filled main levels. Because 
there are no new gaps added on the n*-th step, in this case ~* = is the steady 

state gap-area. 
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a) 



n*(l) 



m*(2) 



m'(l) 



n*(2) 



m*(5) „ ,5, 



0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5 

a 




Figure 3. a) Dependence of to* and n* on a for three r values: r = 1,2,5. Notice 
that relation ( |30| ) holds, b) Dependence of to* and k* on a for r = 5. fe* (diamonds) 
can take on any value from the set {0, 1, 2, .., r — 1}. In all cases a ~ 10 -5 . 



When m* = n* + 1, the last main level is partially filled, and in this case 



A £ = A a ' (on the (n* + l)-st step no new gaps are created) is the steady-state 
gap- area: 



AV = A 



(n*-l) 



(n*) (n*) 

ml 'wl 



(31) 



fc=i 



where k* is defined as the sub-level below which no gaps exist within the main level n*, 



i.e., w£ ^ > for k < k* , and w[" ^ < for k > k* . This results in: 



k* 



where 



[y] + 1 , if y is non-integer 
y , if y is an integer , 



i / l-2a l-a T a"*"' 
V 2ct a r l-a"* T 



(32) 



(33) 



For the m* = n* case we define k* to be zero. Figure |3] shows the values of m*, n* and 
k* as function of a for a few particular r values. 



^.5. Steady-state areas 

Next we analyze the steady-state expression of the area occupied by material B. 
When m* = n* } the steady-state value of the S-area is given by A^ , and for 



m 
A 



n* + 1 by A^' , the complementary of ( pi]) in the unit square, i.e., Ag 



(n*-l) 



(n*) 



- Sfcli m i™ ^ The expression of A ( g L> is computed from (|22|), (|23|), and 



i(n*-l) 



There are two cases in Eq. ( p3| ) corresponding to the generic situation when m is 
non-integer and to the rather special case when u is an integer. They can be treated on 
the same footing by introducing the 'remainder': 
_ j {u} , if u is non-integer , 



if u is an integer 



(34) 



The Active Baker Map 12 
where {x} denotes the fractional part of x (0 < {x} < 1). Thus n* — 1 = u — r. 




Figure 4. The values m* , n* and fc* at r = 3 and er = 10 5 a), and the corresponding 
steady state area A^ ^ of material B, together with A^ , b). Note that the two 
areas coincide for k* = 0, that is when m* = n*. 

From Eq. it follows that a UT = (1 + f/a)" 1 , where v is given in (|2"T|). Using 
this expression, and (p2|), one obtains: 



where 



B V 2(to + (1 - a-) (1 - 2a) J I- (2a) T ' V ; 

1 — a T 1 — a r 

£) = (2a)" rr a + 2 _rT (l-2a) . (37) 

Unfortunately the coefficients C and .D cannot be simplified further, since they contain, 
via r, the fractional-part function which usually hinders analytical manipulations. 
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All these calculations refer to the area of B occupied after n* — 1 active-steps. When 
m* = n* (see Figures |]) and f|a)) Ag — Ag is the steady state area of B. When 
m* = n* + 1 (see Figures [|) and |4]a)) we have to substract from A^ the 'remainder' 
YX=i m( k ^i" ^ This leads to a similar expression as (p5|), but C and D are changed. 
Thus for the steady state area we have 

An*) _ f ^a C + (1 - (2«r) D' _ \ 2a 

~{ 2aa + {l-a-){l-2a) ) 1 - (2a)- ' [ ' 

where C and D' are shorthand notations for 

C = 2 k *~ rr 1 ~ (2a)T (l - a^ 1 -^) + (2a) (1 " r)r , (39) 



1 - a T 

D'= (2a) k *~ rT a + 2 fc *" rr (l-2a) . (40) 



From these we recover (36) and (|37|) with k* = 0, that is, for m* = n*, see figure 4. 



This means that the steady states are described by fl38D , (|39f) , and (^0| ) in all cases. 
Note that the steady state area is a continuous function of a while 1 is not 
necessarily. However, the derivative of A^ with respect to a is discontinuous at all the 
points where k* is discontinuous. 

In this section we have derived the dynamics of the areas covered by materials A 
and B. Our investigations have been restricted, however, to the time instants right after 
the reaction events. In [Appendix A| we include the whole area dynamics of the active 



baker map, i.e., we examine there what happens between subsequent reaction events. 

In the next section we make a connection with the existing heuristical theory of 
active advection and identify the geometric factor for the active baker map. 

3. Comparison with the heuristical theory: derivation of the geometric 
factor 

The heuristical theory presented in the Introduction and also in detail in Ref. @, [| 
derives a recursion relation for the area covered by material B in arbitrary open flows 
just after the (n + l)-st reaction as a function of the area occupied by B just after the 
n-th reaction. This theory is valid for purely hyperbolic systems. In this section we 
show that in the limit of sufficiently small a the phenomenological expressions derived 
in || H are recovered. 

The heuristic formulas were derived with the basic assumption that the material B 
in the steady-state is covering the unstable manifold of a chaotic set. The interesting 
and novel reaction-kinetics is observed in the parameter region (of a, r, Dq and A) where 
the coverage is such that the hierarchical structure of the (fractal) unstable manifold 
is preserved in the system over several lengthscales. For example, if a is too large 
then the fractal filamental structure of the unstable manifold is simply washed out 
during reactions, and thus the reaction-kinetics will be of the common classical type 
of surface reactions [F2]. The effects of fractility of the manifold is seen when a is 



small enough so that the coverage of the filaments of the unstable manifold preserves 
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the filamental structures, in forms of stripes of average width e* < 1. For the baker 
map this means that the number of filaments in the steady state, 2 n * T , is much greater 
than unity: 2 n * T 3> 1. From the calculations above follows that this regime is reached 
for 'intermediate' a values at fixed r and er <C 1. (If a is close enough to either 
or 0.5, n* is close to unity). However, the smaller a is the wider the 'intermediate' 
a-range is where fractal behaviour is important, see the expression of u, Eq. fl2~4]). In 
particular for a = 10~ 5 the number of resolved levels is maximum (n* = 4) in the 
interval a G (0.29,0.49), see figure 4a). In the following we will assume that we are in 
this scaling regime a < 1. 




Figure 5. Comparing the approximate expression ( ]5| ) with the exact (p8|). The solid 
line is Eq. (J3|), the dashed line is Eq. (]4l|). We have also plotted an upper (using g u , 
dotted line) and a lower bound (using ge, dashed-dotted line) for Eq. (^lj). Herer = 3 
and a = 1CT 5 . 



First observe, that in the fraction multiplying 2 U in Eq. ( |38| ) one can neglect <r, since 
a <C 1. Similarly, it is easy to show that 2" ~ (v / ' a )- ln2 / lna (according to (P?|), since 
v ^> 1). It also follows that 2 U 3> 1 and the 1 can be neglected in the curly brackets of 
(|38|). Thus we end up with 

2aD' aq \ 2 -°" 



^■13 



(1 



2a) 



°9 

,kt/(2-D ) 



with D and k given by Eqs. ([12]) and (|13l) , and D' by ([TOD. Here 

±u o (2 J D / )^ 5 o 



2a 



9 



2a 



a 



2a 



(42) 



is a factor which mainly depends on the characteristics of the chaotic baker dynamics. 
Equation ([41]) gives the expression of the steady state 5-area just after a reaction. This 
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coincides precisely with the heuristically derived Eq. (|6]) (see also Eq. (8) of Ref. ||), 
after employing the relation k/(2 — D ) = A. However, while in and in |3|] g appears as 
a phenomenological constant, in the present case of the active baker map, its expression 
can be given explicitely. 




Figure 6. Comparing the geometrical coefficients G, Gi and G u . Here r = 3 and 
a = 1(T 5 . 

Figure |5| shows (on a linear-log scale) for comparison the a-dependence of the steady 
state areas of B for the exact expression of ^ from (^SJ) and the approximation in 
(p|). Note that the two are virtually identical in the scaling regime of a. 

The explicit expression for g is rather complex due to its dependence on the integer- 



part function via D', see Eq. (p}|). It can however be shown that D' is a bounded 
variable. 

First, D' < 1 — a. For the geometric factor this leads to the solely a-dependent 
upper bound, g u : 



2a (\ 



^<^r =9u (43) 

1 — la \ a ) 

Second, a lower bound can be derived for g, as follows. Let F(x) = a{2a)~ x + — 2a)2~ x , 
then D' = F(rr — k*). It can be shown that < rr — k* < 1, and that the function F 
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has a single minimum in the x G [0, 1) interval, at: 



2a / 

Thus, D' > F(xe) = .D^ and therefore g > where ge is given by (|42|) with .D' replaced 

Note that for a — > 1/2, (7 diverges since it is inversely proportional to 1 — 2a. This 
divergence, is however cancelled in the expression for the area by the exponent 2 — Dq. 
In other words, the quantity g 2 ~ D ° is well behaved (1 < g 2 ~ D ° < 4). Let us introduce 
therefore the quantity G = g 2 ~ D °. In figure 6 we show G, Ge and G u as functions of a. 
Figure 5 shows the corresponding area dependence created with both the upper bound 
(dotted line) and the lower bound (dashed-dotted line) . 

Both bounds are independent on r and o~. The two bounds give a rather close 
approximation to the exact area expression, as shown in figure 5. This implies that 
the factor g is weakly dependent on the reaction range a and on the reaction time 
r, corroborating our claim that it is only geometrical, i.e., it mainly depends on the 
parameters of the mixing dynamics (in this case, on a). 

It is also worth comparing the average width of the stripes of material B in the 
baker map to that obtained in the Intoduction [see Eq. (H)] and in Ref. @. The average 
stripe width in the steady state can be obtained as e* = A { b ] l /2( ft *- 1 ) T+fc ^ because there 
are 2 <n *~ 1 ) T+fc * stripes. In the small ff < 1 range this gives 

* „ 2" r 2aD> _ g'a 

£ - 2 (™*- 1 )-+ fc * (1 -a T )(l-2a) 1 - a- ' 1 ' 

where g' is given by 

g' = 2+ Z -?——. (46) 
1 — 2a 

Equation (151) is exactly the same as (|2|) and the steady state of Eq. (6) of Ref. M 
calculated after the reactions. Note that geometrical factor ( [46]) differs from ( fi2]) . The 
reason is that in the heuristic theory of the Introduction and of Ref. || the fractal 
scaling relation Ab ~ e 2 ~ D ° was extended for the steady state stripe width, that is, 
for e*. For this value of e, we are at the edge of the validity of the above mentioned 
scaling, since a crossover to a two-dimensional behavior sets in. This leads to different 
geometrical factors for the areas and for the stripe widths. 



4. Another form of the baker map 

In this section we slightly modify the baker map. One reason to do so is to check whether 
the results obtained so far are robust enough. We show that only the geometrical factor 
is altered in this case, the main results obtained in the Introduction and in Refs. @, |3[] are 
the same. The other reason why another form of the baker map is introduced is based 
on the observation that during the reactions there is material hanging over the vertical 
edges (the edges parallel to the background flow) of the unit square. The behaviour 



The Active Baker Map 



17 



of this material is not described by the baker map in the form of fliCf). Moreover, 
the cutting of these stripes hanging over the unit square is quite unnatural in terms of 
hydrodynamics: there is enough space for the growing of material across the background 
flow, the growth is usually not limited by such a sharp obstacle as the edge of the unit 
square. 

The modified baker map differs from (TO) in the location of the fixed points around 
which the compression takes place. It consists of two steps: 1) the lower and upper 
half of the unit square is compressed in the x direction by a factor of a < 1/2 around 
the fixed point in (1/2 — d/2(l — a), 0) and (1/2 + d/2(l — a), 1), respectively, where d 
is a parameter chosen appropriately, and then 2) the compressed stripes are stretched 
along the y direction (around the same fixed points) by a factor of b — 1/a > 2. In 
other words, after one baker-step we get two vertical (parallel to the y axis) stripes 
centered around (1 — d)/2 and (1 + d)/2. Their centers are at a distance of d. It can 
be chosen such that the stripes do not reach the vertical edges of the unit square, thus 
there remains enough space for the reactions. An appropriate choice is, e.g., d — 1/2. 

The price for this modification is a more involved algebra. In this case, after each 
baker-step there are empty areas (i.e., covered by material A) inside the stripes of width 
a, into which the unit square is mapped. This means that the gaps are wider than in the 
former derivation, and their calculation is not trivial. In fact, it is more convenient to 
follow the distance of the centers of the neighbouring stripes. In a baker-step, all existing 
central distances are multiplied by a, and a new level of the hierarchy is introduced. This 
hierarchy changes only if there are gaps disappearing. Note that the stripe dynamics of 
the areas covered by material B is much simpler than in the former baker model: all 
stripes of material B have equal widths. 

Without going into details, we summarize the results obtained with this modified 
baker model for the steady state areas. First let us introduce 

d(l-2a)(l-a T ) 

a C R = > . (47) 

This is a critical value of the reaction distance for any choice of a and d; if a > <7cr 
then there is material B hanging over the vertical edges of the unit square in the steady 
state. For a < <tcr, which we assume in what follows, we obtain a steady state where 
a Cantor-like hierarchy is produced with N different levels, where N is given by 

N = f [z] + 1 , if z is non-integer , 
1 z , if z is an integer , 



where 



In 



Hi) ' 

The steady state stripe width of material B after the reactions is 



(49) 



— T + ¥±. (50) 
— a T 1 — a 
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This is to be compared to the steady state stripe width (0), using (|12|) and (|TT|), 
which hold for the modified baker map as well. From this equation and from ( |50D we 
obtain 

/ l — a T a N d 

9 = 2 + -. • 51 

1 — a a 

In case of this modified baker map we have also obtained a correction to the g — 2 value, 
which is due to the overlapping of stripes during the reactions. The geometrical factor 
depends mainly on parameter a of the baker map, as from (|47|) and (^) we obtain the 
bounds 

2a . 2 

2 + —a^^ 2 + —a ■ < 52 > 
The total area covered by material B is 

A B = 2»e> = 2»(-^ + -^) . (53) 
\ 1 — a T I — a) 

This is, again, in coincidence with that obtained in Eq. (8) of 0, but the geometrical 
factor is different from ([51]), it is 



1 - a T 2 N+T+1 aa T (2a) N+T d\^ 

9 = — — ; — + V • ( 54 ) 

aa T \ 1 — a T 1 — ay 

The deviation between ~g and g' is again the result of the extension of the fractal scaling 
to the range of e*, to the steady state stripe width. 




In Figure we compare the four geometrical factors d42|) (solid line) and (|54|) (dotted 
line), and (|46D (dashed line), (|5l| ) (dashed-dotted line) as function of a. Although the 
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average trend of the curves is the same, there are large deviations due to the fact that 
g' and ~g' are discontinuous functions. The deviations, however are less pronounced in 
the scaling region where the fractal structures of the maps are better resolved. This 
leads us to the conclusion that the geometric factor strongly depends on the particular 
realizations and relative positions of the filaments even when these realizations share 
the same fractal and dynamical properties (D Q , A). 



5. Discussions and outlook 



In this paper we have investigated the effects of a fractal geometric catalyst supplied 
e.g. by the Langrangian dynamics of an underlying flow on the reaction kinetics of 
autocatalytic reactants transported in this flow. The reactive particles were assumed to 
simply follow the local velocity of the fluid element at their position and act as simple 
passive tracers that do not modify the properties of the flow. For mixing dynamics we 
have considered the simplest possible dynamics, namely two variants of the baker map, 
in order to be able to generate an analytic description of the reaction process. Even 
though the baker map is the simplest map with chaotic mixing, when reactions are 
superimposed on the transport dynamics, the analytical calculations become involved 
since the reaction range a introduces a finite nonscaling cutoff which breaks the strict 
selfsimilarity. In the limit of slow reactions one is able to recover from these exact 
formulas the heuristic expressions derived for general flows, thus validating the latter 
approach, and showing that indeed, the fractal backbone does enhance the productivity 
of reactions, and acts fractal catalyst. 

A crucial point in the consistency between the heuristic and exact approach is played 
by the so-called geometrical factor g (or its variants g',g,g'). By definition, it should be 
a parameter reflecting mainly the geometrical properties of the mixing dynamics (of the 
hydrodynamical flow). In the heuristic theory it is postulated to be independent of the 
chemistry parameters a and r. The consistency is achieved by showing that it is only 
weakly depending on these parameters, in the scaling limit at least. This property follows 
from the close bounds (|43D and (|52|) on the geometrical factor(s) depending on a only. A 
direct plot of the geometrical factors as a function of the reaction range a (figure 8) also 
supports this property and suggests a weak, approximately periodic dependence on the 
logarithm of the reaction range, over several orders of magnitude. Note, that the actual 
value of the geometrical factor strongly depends on the two independent global chaos 
parameters (A, Do) which determine the chemical kinematics (see Introduction), but 
might also depend on other aspects of the chaotic dynamics (as e.g., on the parameter 
d of the modified baker map). 

Various improvements and extensions to the present theory are possible: 1) one 
can consider [TB|, |I3J the more realistic situation of non- simultaneous reactions, i.e., 
the particles have a random phase dispersed among them. We expect that this only 
leads to a more rough, fluctuating A — B interface along the manifolds, however these 
fluctuations are bounded due to the fact that only correlations of finite length can 
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Figure 8. The dependence of the geometric factors g, g, g' , g' on the reaction range, 
a. Here a = 0.4, r = 3, and d = 1/2. 



develop along the interface during the finite time lag r, and thus it will renormalize 
the front velocity; 2) the interacting particles (which can be for example phytoplankton 
[14|) may locally modify the advection dynamics due to inertia and finite size effects. 
In this case the phase space of the problem increases in dimensionality by the number 
of velocity components of the particles since the latter no longer is the same as the fluid 
elements' velocity; 3) and at last but not the least, various other chemical reactions 
should be monitored, besides the autocatalytic ones presented here. 

Finally, it is worth connecting the baker map results with the general description 
given in the Introduction. To this order, we mention that in a map of linear extension 
L, the reaction range is oL. Analogously, the dimensional reaction time is tT, where T 
is the temporal period over which the baker map is acting. The reaction velocity v r is 
then 

oL 

v r = — . (55) 

Taking the continuum limit of the reaction (see Introduction) means that both the 
reaction range and time are short on the corresponding macroscopic scales. Since L is 
interpreted as a hydrodynamical length, this implies tr < 1 as used in the baker map. 
On the other hand, since r was taken to be of order unity, tT can only be short if the 
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duration of the map is much shorter than the characteristic hydrodynamical time. Even 
without further specifying the choice of T, we can derive a useful condition in terms of 
a dimensionless number. 

In the case of reactions in fluid flows with chaotic mixing properties we have three 
distinct levels of time-continuous dynamics built on top of each other: 

• the Eulerian description of the hydrodynamics, characterized by one main 
parameter, the Reynolds number 

Re = — , (56) 
v 

where U is a typical velocity, L is the typical length scale and v is the kinematic 
viscosity; 

• the Langrangian description of the dynamics, characterized by the fractal dimension 
of the 'transport manifolds' (unstable manifolds), and by the contraction dynamics 
on them, 

D = 2 - j , A = \(Re) , k = n(Re) ; (57) 
A 

• the reaction kinetics, described by a 'reaction' Peclet number, 

(Per)- 1 = ^- , Pe r = Pe r (X(Re)). (58) 
XLi 

Reaction equation ([/]) can be made dimensionless if the area is measured in units 
of L 2 , and the time in units of the characteristic lifetime of (the Lagrangian) chaos, 
1/k=1/(A(2 -£>„)): 

A B = -A B + g ^ (A B y P ■ (59) 

The active baker map's Peclet number is 

Here we have used that the dimensional Lyapunov exponent of the map is — In a/T. 
With our typical choice of a = 1CT 5 , r = 1, 2, 3, and a on the order of 1/4, the reaction 
Peclet number is much larger then unity. 

It is worth giving the definition of the traditional Peclet number [16] which measures 
the importance of diffusion relative to the flow: 

(Pe)- - £ - a (6!) 

where — D/L is the charactersitic diffusion velocity. Large Pe implies weak diffusion. 
As suggested by the analogy between the traditional and the reaction Peclet number, 
the reaction range plays a similar role as a diffusion length and the reaction velocity 
v r as D/L. Furthermore, if molecular diffusion is present but not very strong, it only 
renormalizes the reaction front velocity and all the formulas shown here are valid with 
v r replaced by an effective reaction front velocity M. 
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Note that the direct reaction analog of the Peclet number would be v r /U. It is 
worth pointing out however that instead we have 

(Pe ) _1 - — — - — L y a P unov time / 62 ) 
U XL U hydrodynamical time 

which is rather a Lagrangian dynamics based quantity. The difference between Pe r 
and its hydrodynamical analog, v r /U, might be pronounced if the hydrodynamical and 
Lyapunov times are not on the same order of magnitude. 

The limit of slow reactions implies v r <C XL, i.e., Pe r ^> 1. It is in this limit 
where pronounced fractal product distributions are expected. The steady state product 
area A* B given by Eq. (|9|) is by a factor (Pe r ) Do_1 ^> 1 larger than it would be in 
a nonchaotic flow with Dq = 1. Thus we find that a large value of the reaction 
Peclet number introduced by Eq.flSSD provides a general, neccessary condition for the 
existence of filamantal product distributions observable over several orders of magnitude 
in reactions taking place in chaotic flows. 
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Appendix A. Area-dynamics in the active baker map 

In Sec. we described the hierarchy of the stripes formed in the active baker map just 
after each reaction, calculated the total area occupied by the materials A and B as a 
function of the number of active- steps, and proved that a steady state sets in after a given 
number of active-steps, m*. In this appendix we present the full temporal dynamics of 
the areas occupied by A and B (not only as a function of the active-steps). 

In order to account for the 'intermediate' steps as well, we will extend a little bit 
the notation introduced in Sec. |[ First we have to observe that once a reaction has 
been completed and new passive baker-steps are taken, the existing hierarchy at the end 
of the last reaction is doubled and mapped on a level one step lower in the hierarchy 
after a passive baker-step, while at the top of the hierarchy there are reaction-free gaps 
appearing. If p denotes the number of reaction-free baker-steps taken just after the last 
reaction, the top gap of width 1 — 2a appears at p — 1, the additional two gaps of width 
a(l — 2a) appear at p = 2, etc., until at p = r the reaction-free tree becomes completed. 
All the rest of the gaps coming from earlier reactions are lying under this newly created 
'baker tree' of length r. Next a reaction occurs and the whole process repeats itself. 
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a=0.45,T=l^ 
a=0.45, x=2 --■*- 




5 10 15 20 25 

t=nx+p 



Figure Al. The evolution of the area of material B shown on a log- linear scale for 
three values of r and two values of a. In all cases a = 10~ 5 . 



Let us denote the widths of the gaps coming from p baker-steps by w\ , £ G 
{1,2, ..,p}. We have wf'^ = a^ _1 (l — 2a) and the corresponding multiplicities are 
r7 4°' P ^ = 2^ _1 . We shall denote the widths and multiplicities of the gaps resulted 
from previous n reactions and the additional p passive baker-steps by w^' p \ and mjjf , 
respectively, where k = 1, 2, .., r and j = 1, 2, .., n (the sub-trees coming from previous 
reactions are all completed). To connect with the notation from Sec. observe 
that w^' ^ = w\. . We have = a p w^\ and the corresponding multiplicity of 

mi' p) = 2Pm£\ k = 1, 2, .., r and j = 1, 2, .., n. 

When the (n + l)-st reaction occurs we obtain the widths where k = 1, 2, .., r 
and j = 1,2, ..,n + 1, with tojjf = itf* — 2cr, mj^ = m^'^ (and it also leads to 
Eq. fl2~0|)). The area of material A (B) after n reactions and p reaction-free baker-steps 
will be denoted by A%' p) (Ap p) ). By definition, A^ 0) = A^ (Ap 0) = A^). The 
areas can be calculated as follows: 

A%' p) = J2 mf' p) wi°' p) + £ £ m ( t p) w^ p) = 1 - (2a)' + (2a)" (A.l) 

i=l j=l k=l 

where A^ is given by (^). The above expression is valid in the regime where no gap 
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fillings have occured yet. In terms of the 5-material, using flUf), ( |A. 1| ) is equivalent to 

Ap p) = {2a) p A^ = e- K 04 n) (A.2) 

just as expected, since between two reactions there is an exponentially emptying 
dynamics. The total time is t = nr + p, p = 1,2, ..,r. When p = r, t = {n + l)r 
and a reaction occurs instantaneously. 

Figure |A1| shows the t-dependence of the area covered by the material B for three 
r values, for given o [o = 10~ 5 ) on a log-linear scale at two a values, a = 0.45 and 
a = 0.32. Notice the fact that the actual time evolution of the areas follow closely the 
simple exponential emptying dynamics. Certainly, between two consecutive reactions, 
the emptying is exponential, the deviation from the exponential in the overall run of 
the curves occurs in jumps at reaction events. 

In order to observe a longer temporal evolution before stationarity sets in, one needs 
to choose a very small a. On the other hand, this means that the overall fattening during 
a reaction is very weak, and it is in fact a 'perturbation' at early times. It becomes 
'visible' just before the steady state. 
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